Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN WT^ style file vl.4) 



Non-linear evolution of the tidal angular momentum of 
protostructures II: non-Gaussian initial conditions 

Paolo Catelan^^ and Tom Theuns ^ 

^ Department of Physics, Astrophysics, University of Oxford, Keble Road, Oxford OXl 3RH, UK 
^ Theoretical Astrophysics Center, Juliane Maries Vej 30, DK-2100 Copenhagen, Denmark 

1 February 2008 



^ . ABSTRACT 

^ il The formalism that describes the non-Unear growth of the angular momentum L of 

■ protostructures from tidal torques in a Friedmann Universe, as developed in a previous 

paper, is extended to include non-Gaussian initial conditions. We restrict our analysis 
here to a particular class of non-Gaussian primordial distributions, namely multiplica- 
tive models. In such models, strongly correlated phases are produced by obtaining the 
gravitational potential via a nonlinear local transformation of an underlying Gaussian 

' random field. The dynamical evolution of the system is followed by describing the 

trajectories of fluid particles using second-order Lagrangian perturbation theory. In 
the Einstein-de Sitter universe, the lowest-order perturbative correction to the vari- 
. ance of the linear angular momentum of collapsing structures grows as t^^^ for generic 

' non-Gaussian statistics, which contrasts with the growth rate characteristic of 

■ Gaussian statistics. This is a consequence of the fact that the lowest-order perturba- 
I five spin contribution in the non-Gaussian case arises from the third moment of the 
. gravitational potential, which is identically zero for a Gaussian field. Evaluating these 

Oh' corrections at the maximum expansion time of the collapsing structure, we find that 

Q . these non-Gaussian and non-linear terms can be as high as the linear estimate, without 

^ ' the degree of non-Gaussianity as quantified by skewness and kurtosis of the density 

. field being unacceptably large. The results suggest that higher-order terms in the per- 

turbative expansion may contribute significantly to galactic spin which contrasts with 
the straightforward Gaussian case. 

' Key words: galaxies: formation - large-scale structure of Universe 

1 INTRODUCTION 

It is often assumed that protostructures formed as primordial density fluctuations are amplified by gravity. The initial density 
field is generally assumed to be a Gaussian random field and the amplitude of fluctuations on all scales is then characterised 
by the density power spectrum. An appealing consequence of this gravitational instability model is that patches of matter are 
spun- up, as density inhomogeneities exert tidal torques on collapsing objects (Hoyle 1949), thereby providing an explanation 
for the spin of galaxies. The amount of angular momentum tidally induced on a given patch of matter depends mainly on its 
shape and scale R (mass A/), and to a lesser extent also on its environment. Simple arguments based on linear gravitational 
instability theory lead to values for the angular momentum of the Milky Way which compare favourable to observational 
estimates (Peebles 1969). The statistical distribution of the spin of collapsing objects has been obtained by studying tidal 
torques on (high) peaks of the linear Gaussian density field (Hoffman 1988; Heavens & Peacock 1988), in the formalism of 
Peacock and Heavens (1985) and Bardeen et al. (1986; BBKS). More recently, we adopted the Lagrangian approach suggested 
by White (1984) and the BBKS peaks formalism to derive the probability distribution of the rms spin L (Catelan & Theuns 
1996a): we showed that the distribution of this variable is practically indistinguishable from that of the modulus jL| as 
obtained by Heavens & Peacock (1988), indicating that the rms spin is well suited to characterise the statistics of rotational 
properties. In addition, we used the same formalism to derive the distribution of the specific angular momentum L/M versus 
M on galactic scales (Catelan & Theuns 1996a), which allows a direct comparison with the observational data reported in 
Fall (1983) as well as with the numerical simulations of Navarro, Frenk & White (1995). 
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The above investigations have been restricted to the linear regime, during which the galaxy spin grows proportional to the 
cosmic time t (Doroshkevich 1970; White 1984). We have carried out a theoretical analysis of how the tidal angular momentum 
evolves during the mildly non-linear regime (Catelan & Theuns 1996b; Paper I). The formalism employs pcrturbative solutions 
of the Lagrangian dynamical equations describing the motion of the fluid, and assumes Gaussian initial conditions. The 
Lagrangian approach shows to be very well suited to treat the non-linear evolution of the galactic spin, because it is powerful 
in describing the growth of the mass-density fluctuations on the one hand (Zcl'dovich 1970a, b; Buchert 1992; Bouchct ct 
al. 1992; Catelan 1995); and on the other hand, the usual difficulty of inverting the mapping from initial Lagrangian coordinates 
q to Eulerian coordinates x is completely by-passed. The latter is because the angular momentum L is invariant with respect 
to an Eulerian or Lagrangian description. In spite of the fact that the leading spin corrections grow faster (oc t^^^'^) than the 
linear term, the main result of the perturbative investigation as reported in Paper I is that the predictions of linear theory are 
rather accurate in quantifying the evolution of the angular momentum of protostructures before collapse sets in: non-linear 
corrections increase the rms spin by a factor ~ 1.3 for the CDM power spectrum on galactic scales. 

In this paper we extend the perturbative investigation of Paper I by considering non-Gaussian initial conditions. Specifi- 
cally, we want to address the question whether allowing the primordial matter field to deviate from the Gaussian distribution 
modifies the subsequent evolution of the angular momentum during the mildly non-linear regime. 

A first problem we have to face is the wide variety of possible non-Gaussian probability distributions. Primordial non- 
Gaussian fluctuations are predicted in many different cosmological scenarios. Topological defects remaining after an early 
phase-transition (Kibble 1976) such as cosmic strings (Zel'dovich 1980; Vilenkin 1981; Turok 1984; Scherrer, Melott & 
Bertschinger 1989), monopoles (Bennet & Rhie 1990) or global textures (Davis 1987; Turok 1989; Turok & Spergel 1990) are 
examples of models whose statistics are not Gaussian. Modifications of the inflationary scenario leading to phase correlations 
on cosmologically relevant scales have been discussed (Allen, Grinstein & Wise 1987; Kofman & Pogosyan 1988; Salopek 
& Bond 1991; Salopek 1992): these are non-Gaussian as well. Following suggestions of Kofman et al. (1989), Moscardini et 
al. (1991) and Weinberg & Cole (1992) considered non-Gaussian models obtained by performing non-linear transformations 
of an underlying Gaussian random field. In particular, this last class of non-Gaussian models is interesting in that one is able, 
by tuning free model parameters, to reproduce some of the particular properties of, e.g., the texture-seeded Cold Dark Matter 
model (Gooding, Spergel & Turok 1991; Park, Spergel & Turok 1991) or the cosmic explosions scenario (Ostriker & Cowie 
1981). 

Non-Gaussian perturbations constitute a more general statistical model than Gaussian ones and can be adopted to com- 
pute cosmological observables (see Wise 1988), such as spatial galaxy correlation functions (Matarrese, Lucchin & Bonometto 
1986; Scherrer & Bertschinger 1991), peculiar velocity correlation functions (Scherrer 1992; Catelan & Scherrer 1995; Moessner 
1995), expected size and frequency of high density regions (Catelan, Lucchin & Matarrese 1988; Matsubara 1995), hotspots 
and coldspots in the cosmic microwave background radiation (Coles & Barrow 1987; Kung 1993), and higher order temperature 
correlation functions (Gangui et al. 1994; Moessner, Perivolaropoulos & Brandcnberger 1994). 

We limit our analysis in this paper to multiplicative non-Gaussian processes, which have been analysed recently with 
AT-body simulations in the context of a Cold Dark Matter cosmology (Messina et al. 1990; Moscardini et al. 1991; Matarrese et 
al. 1991; Weinberg & Cole 1992). In such models, the gravitational potential is obtained via a local non-linear transformation 
of an underlying Gaussian random field. These models are interesting because the procedure provides strongly correlated 
phases, while at the same time preserving the standard form of the initial power spectrum. Consequently, Gaussian and 
non-Gaussian models with identical power spectra can be compared to draw conclusions on the effects of departing from the 
more usual assumption of a normally distributed primordial random field. Infiatiori-generated non-Gaussian fiuctuations are 
expected to be of the multiplicative type (Matarrese, Ortolan & Lucchin 1989; Kofman et al. 1990; Barrow & Coles 1990). 
An important example of a multiplicative process is the lognormal statistic, which has recently been proposed as a simple 
phcnomcnological model for describing the (present day) non-linear density field (Coles & Jones 1991; see also Bernardeau & 
Kofman 1995). 

In Section 2 we review how Lagrangian perturbation theory can be used to obtain corrections to the linear angular mo- 
mentum due to non-Gaussian initial conditions. In Section 3 we introduce the non-Gaussian statistics that we will investigate. 
In the next Section, we combine the obtained results to compute the spin corrections for the different non-Gaussian statistics. 
Finally, Section 5 summarises our findings. Technical Appendices contain details on different aspects of the computations. 



2 NON-LINEAR EVOLUTION OF THE TIDAL ANGULAR MOMENTUM 

Let us assume that the behaviour of matter on scales smaller than the horizon is similar as that of a Newtonian pressureless 
and irrotational self-gravitating fiuid embedded in an expanding universe with arbitrary density parameter fl. For the sake 
of simplicity, we consider the case of zero cosmological constant and leave to the interested reader the corresponding general- 
isation. Additionally, luminous objects like galaxies and clusters are assumed to grow through gravitational amplification of 
primordial positive density fluctuations 6 of the coUisionless fluid. 



© 0000 RAS, MNRAS 000, 000-000 



Evolution of the tidal spin from non- Gaussian initial conditions 3 



In what follows, comoving coordinates are denoted by x and physical distances by r = a(t)x, where a{t) is the scale 
factor and t the standard cosmic time. The dynamical fluid equations in a non-fiat Friedmann universe simplify considerably 
if, instead of t, one adopts the temporal coordinate r defined by dr = a"^ dt (Shandarin 1980). In fact, in terms of r, 
the peculiar velocity and the peculiar acceleration are given by dx/dr = x = a(r)u(x, r), and d^x/dr^ = x = g(x,r). 
As discussed in Shandarin (1980), this dimensionless time r is negative and the initial cosmological singularity at t — 
corresponds to r = —oo. Furthermore, the infinity of the cosmic time, t = -|-cx3, in the open models corresponds to r = — 1, 
whereas t — +oo in the critical Einstein-de Sitter universe corresponds to t = 0. Finally, the contraction phase in the closed 
models starts at r = 0. In terms of the density parameter 57, one has r = —\/—k (I — f2)^^''^, where k is the curvature 
constant (A; = — 1 for open universes and k = \ for closed universes). The case f2 = 1 (fc = 0) is a singular point for the latter 
transformation and in this case we take r = —(?>i)~^^^, which corresponds to using o(t) = (3t)^'^^ or to = t/a^^'^ = 1/3, so 
defining the unit of time. The scale factor a(r) may then be written for all Friedmann models as a(r) = (r^ + k)~^ . 

The linear evolution of angular momentum of proto-objects is most easily analysed using the Zel'dovich (1970a, b) 
formulation (see White 1984; Catelan & Theuns 1996a) and the mildly non-linear spin growth is most easily analysed using 
Lagrangian perturbation theory. We recall that the Zel'dovich approximation coincides with the linear Lagrangian description. 

We will essentially adopt the formulation of the Lagrangian gravitational theory for a coUisionless Newtonian fluid as 
presented in Catelan (1995; see also references therein) but note that in the present paper, as we did in Paper I as well, the 
variable r has opposite sign and the growth factor of the linear density perturbation is normalised differently. An alternative 
formulation of the Lagrangian theory may be found in Buchert (1992). 



2.1 Basic tool: Lagrangian theory 

In the Lagrangian formulation, the departure at time r of mass elements from their initial position q is described in terms of 
the displacement vector field S, 

x(q,T) = q+S(q,r) . (1) 
The trajectory S(q, r) satisfies the Lagrangian 'irrotationality' condition and the Poisson equation given by (Catelan 1995) 

£(1/37 [(1 + V ■ S) 5/3^ — S/la + S^ct] S^a = , (2) 

[(1 + V ■ S) 5„fi - S^i3 + S^p] Sp^ = a(r)[J(q, r) - 1] , (3) 

respectively, where e^p-y is the totally antisymmetric Levi-Civita tensor of rank three (ei23 = 1), the symbol Sap indicates 
the Kronecker tensor, and summation over repeated Greek indices (where a — 1,2,3) is understood. In these equations, 
a(r) = 6a(r) and J = 1/(1 + 5) is the determinant of the Jacobian of the mapping x ^ q. This determinant J is non-zero 
until the first occurrence of shell-crossing (see, e.g., Shandarin & Zel'dovich 1989). In addition, Sai3 = dSa/dqp, V = Vq, and 
S^0 denotes the cofactor of Sap- In general, Soip is not a symmetric tensor: Sap = Spa if, and only if, the Lagrangian motion 
is longitudinal in which case S is the gradient of a potential. 

The equations (^) and (j^ are the complete set of dynamical equations for the vector field S, which describes the 
trajectory of massive particles of a coUisionless fluid embedded in an arbitrary Friedmann universe. We briefly summarise 
their perturbative solutions (up to second-order in the linear displacement) in the next subsection. 



2.2 Lagrangian perturbative solutions 

The master equations (^ and ^ are of third-order and non-local in the displacement S (see the discussion in Kofman & 
Pogosyan 1995) and it is undoubtedly very difficult to solve them rigorously. Perturbative solutions can be found by expanding 
the trajectory S in a series of which the leading term corresponds to the Zel'dovich approximation. Specifically, S — E„S„, 
where Sn = 0{S") is the n-th order approximation. Here, we will need Si and S2 only. For brevity we use the symbol 
0(r) = ln[(T — l)/(r -I- 1)]^^^ for the open universe case {k — —1) and A(r) = arctang(l/r) — —iO{iT) for the closed universe 
case (k — +1). We neglect decaying modes. 

2.2.1 First-order approximation: 

The first-order solution to equations (^ and (^) is separable in space and time and corresponds to the Zel'dovich approximation 
(Zel'dovich 1970a, b): 

Si(q,r) = D(r)S(^)(q) = i3(r)V^(^'(q) ; (4) 
the function D(t) is the growth factor of linear density perturbations and is given by: 
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( |{l + 3(r2-l)[l+re(r)]} if Q< 1 
D{t) = I if f7 = 1 (5) 

[ f{-l + 3(r2 + l)[l-rA(r)]} if f7 > 1 . 

The solution for the closed models can be obtained from that for the open models by substituting in the latter r by ir and 
reversing the sign to make the growing mode positive. Note that, in contrast to Bouchet et al. (1992) and Catelan (1995), we 
normalised D{t) according to the suggestion of Shandarin (1980): the coefficient 5/2 is such that D{t) — > in the limit 
r — » — oo, which coincides with the Einstein-de Sitter case. D{t) for the different universes is plotted in Fig. 1 of Shandarin 
(1980) and in Fig. Al of Paper I. The function tp^^^ci) is the (initial) gravitational potential. For later use we define its Fourier 
transform, i/;'^'(p) = J dq-ip^^'' (q) e~^^''^ , where p is the comoving Lagrangian wave vector. The Fourier transform of the 
linear density field, (5'-^-'(p,r) = D{t)5i{p), is related to '4'^^\p) via the Poisson equation, ^^^\p) = p^^Si{p) . 



2.2.2 Second-order approximation 

The second-order solution is also separable and describes a longitudinal motion in Lagrangian space: 

S2(q,r) = i5(r) S(^)(q) = iJ(r) W^*^' (q) . (6) 
The second order growth rate E{t) is 

[ -f -^(r^-l){l + re(r) + i[r+(r2-l)e(r)]^} if Q < 1 
£(r) = i -fr-* iff7 = l (7) 

[ -f + ^(r' + l){l-^A(r)-i[r-(r^ + l)A(r)]'} if f7 > 1 . 

The second-order non-fiat solution was derived previously by Bouchet et al. (1992). The extra factor 25/4 of the present 
formulation is due to the different normalisation of the first-order solution D. An excellent approximation of the second-order 
growing mode is _E « —jD^ (see e.g. Fig. 5 of paper I). In the limit r — > — oo one has E — — Ir""* which corresponds to the 
fiat case. The Fourier transform of the second-order potential i/)'-'^^ is (Catelan 1995): 

= I ^|^[(2-)='fc(pi + P2-p)]^<^'(pi,P2)^'^'(pi)^(^'(p2), (8) 
where we have defined the symmetric kernel 

«*^''(Pl,P2) = ^ [p?P2 - (Pl- P2)^] , (9) 

which describes the second order non-linear corrections to the trajectory of the fluid elements. 

The perturbative expansion for S corresponds to a Taylor series in the variable D{r) — in the Einstein-de Sitter 
universe but this is no longer rigorously true in a non-fiat universe. However, since the higher-order growth factors can be 
approximated exceedingly well by powers of D, the expansion in the non-fiat case is still 'close' to a Taylor expansion (see 
the discussion in Paper I). We proceed by considering how corrections to the linear displacement Si translate into corrections 
to the linear angular momentum. 



2.3 Non-linear spin dynamics 

In this section, we briefiy summarise the perturbative approach to the spin evolution. As stressed in Paper I, the angular 
momentum of the matter contained in the volume V{t) at a given time r may be written equivalently as an integral over the 
corresponding initial volume F: 

L(r)=r,o^dq[q+S(q,r)]x^^^. (10) 

Here, rjo = a'^pt, where pi, is the mean background density. This procedure enables us to apply the Lagrangian description of 
Newtonian gravity previously reviewed. The linear regime (Zel'dovich approximation) has been fully analysed in this way by 
Doroshkevich (1970), White (1984) and Catelan & Theuns (1996a), whereas its Eulerian counterpart was studied extensively 
by Heavens and Peacock (1988). We can extend the linear Lagrangian analysis of the evolution of the angular momentum 
L(r) to the non-linear regime by applying perturbation theory to equation (0). Perturbative corrections to S(q, r) (Bouchet 
et al. 1992; Buchert 1994; Catelan 1995 and references therein) then give perturbative corrections to L(r): formally 

oo oc h „ 

L(r) = Yl ^Y.T.^'^r'' ^^■^'^'"^ " ^^^^i^lHll) , (11) 



h=0 j=0 



r 
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with So = q, hence L'"' — 0. To calculate the lowest-order corrections to the ensemble average (L'^), we need to compute 
corrections to L up to second-order. After reviewing briefly the results of the linear theory, we summarise the final expression 
of the correction L'^-*. 



2.3.1 Linear approximation 

The linear Lagrangian theory corresponds to the Zel'dovich approximation and the first-order term in equation ( pT] ) is given 
by: 

L^^\r) = 770 D{r) dqq X V^^^\ci) . (f2) 

Assuming that 7/;'-^'(q) can be adequately represented in the volume F by the first three terms of the Taylor series about the 
origin, q = 0, each component L^a\t) may be written in compact form as (White 1984; Catelan & Theuns 1996a): 

dp 

(2^ 



L^i' {t) = b{T) e^p^ O^t' J-i = -D{t) e,^^ J^^ / p^ p^ W{pR) i>^^^ (p) , (13) 



where we introduced the deformation tensor at the origin, O^^J = ©^^J(0) = 9/39^1/'^^' (0), and the inertia tensor of the mass 
contained in the volume F, = rjo dq^o- q-t • In addition, the field tp''^^ is now assumed to be filtered on scale R using the 
smoothing function Wr, whose Fourier transform is W{pR). Equation (^) shows that the linear angular momentum L'-^-' is in 
general non-zero because the principal axes of the inertia tensor Ja/s, which depend only on the (irregular) shape of the volume 
F, are not aligned with the principal axes of the deformation tensor 2?^^^, which depend on the location of neighbour matter 
fiuctuations. The temporal growth of the tidal angular momentum induced by this misalignment is completely contained in 
the function D{t), which behaves as D{t) = — 2r~'^ ~ t in the Einstein-de Sitter universe, as first noted by Doroshkevich 
(1970). Finally, if F is a spherical Lagrangian volume, then L^ ^ ~ ^a/Sj'D^p^ = 0. Consequently, the matter contained initially 
in a spherical volume does not gain any tidal spin during the linear regime (see also the discussion in White 1984). 



2.3.2 Second-order approximation 

The second-order term in equation ( pT| ) involves the second-order displacement S'^': 

L^''W = '70 / dqq X ^ = VoE{t) / dqq X VV''"(q) • (14) 



J r Jr 

Note that, since E oc r~*, one has E oc hence the second-order terms grows oc t^''^ in the Einstein-de Sitter universe. 
This growth rate was first derived by Peebles (1969). If we represent ilj^^\q) in F by the first three terms of a Taylor series, 
as we did before for tp^-^\ we obtain for the a-component (see Paper I): 

L^^\t) = E{r)e^p^vflj^-, 

T?( \ T ( d-pidp2 (Pl + P2)ct (Pl 4- P2)7 ~ , I 12), N 7(1). N 7(1)/ N 

= EirJeap-yJai -TTTTe 1 ; 12 Pi + P2K?) '(pi, P2) i/^^ '(pi) i/-^ ' P2) , (15) 

J (2vr)t> |Pi + P2r 

where O^^^ = I?^'^'(0) — di}da-'ip^'^\0) is the second-order deformation tensor. In addition, the second-order potential -f/)'^' in 
equation (^^ is now explicitely smoothed on scale R by the filter W(pR). The filtering of the field t/;^^-' on an appropriate scale 
refiects the restriction to the mildly non-linear evolution of the proto-object, since the strongly non-linear couplings between 
different modes are filtered out. Note that the non-linear dynamical evolution modifies only the deformation tensor and not 
the inertia tensor. Furthermore, if F is a sphere, then again L'^-* = 0. This is in contrast to Eulerian perturbation theory 
since the angular momentum of an Eulerian sphere does grow in second-order perturbation theory (Peebles 1969; White 1984). 

For a single collapsing region enclosed in a volume F it is enough to evaluate equation ( pi] ) at the time of maximum 
expansion tm to compute its final angular momentum. After tm the angular momentum essentially stops growing since the 
collapsed object is much less sensitive to external tidal couplings (Peebles 1969). However, highly non-linear interactions 
typically occuring after the maximum expansion time may lead to a significant redistribution of angular momentum in the 
final object in a complicated way (see discussion in Catelan & Theuns 1996a and references therein). 

In order to compare the theory against statistical results obtained from A''-body simulations or from observations, it is 
useful to compute the variance of the angular momentum of the object for an ensemble of realisations of the gravitational 
potential random field xjj'^^K This programme is carried out in the next section. 
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3 NON-GAUSSIAN RANDOM FIELDS: SPIN ENSEMBLE AVERAGES 

We simplify the previous results by considering the expectation value over the ensemble of realisations of the non-Gaussian 
random field i/j'^-* of the square of L, (L^). Note that one can relatively easily obtain the probability distribution for the 
linear angular momentum (Heavens & Peacock 1988; Catelan & Theuns 1996a). However, such a calculation would be far 
more complicated for the non-linear contributions, even in the Gaussian case. We will follow the same procedure adopted in 
Paper I, i.e. we neglect the correlations between the inertia tensor and the gravitational field i/;'^'. 
Taking into account the mildly non-linear corrections one has 

(L") = (L*')') + 2(L<') ■ L*')) + 0(t-^) . (16) 

The linear term (L''^-'^) has been discussed extensively by Catelan & Theuns (1996a). Since it involves only the second-order 
moment of the underlying distribution (i.e. the power spectrum), its formal expression is valid for all gravitational potential 
fields, both Gaussian and non-Gaussian. We report its explicit expression below (see Eq. 49). Furthermore, it is easy to 
understand from equations (|l^) and ( p^ ) that (L^^' • L'^') is the lowest-order perturbative correction since it corresponds to 
an integral over the bispectrum B^^^ of the initial gravitational potential field i/)'^': 

(27r)^ Sn{pi +P2 + P3) B[^\pi,p2) = (^''^(pi) ^'''(pa) ^^'''(pa) ) . (17) 
Here, Sd is the Dirac delta function. Specifically, one gets 

{L(^) • L(2)) = D{r) E{r) e^^^ e,^,y J.^ J,,.,, { , (18) 
where 

(vf^vfl) = -/ ^ ^'^^ -^^^^^ ^''^ ^ Ip! + P2r'^'' ^'^^ ^'^'^^ ^'''^^"^^^ [W(|pr + P.|i?)]^<(p.P2) 

= (19) 

(2) 

which defines the kernel /C^^^,^, describing the effects of the 2-order dynamics. To maintain a compact notation, we have 
introduced the operation "©" defined for integrable functions J- and Q as 

^ee = y"^g^^(pi,P2)e(pi,P2) . (20) 

Saturating the tensors' indices in equation (|l^ ) and using the isotropy of the universe, it is possible to write equation (|l^) in 
full generality as (see Appendix A for details) 

(L(i) . l(^)) = A D(r)Eir) (/i? - 8^2) S(i?) , (21) 
where 

E(i?) = 15 (oiV Vg^ ) = 15 IC^^ly (BB^^K (22) 

Here, and fi2 are the first and the second invariants of the inertia tensor J (section 3.2 below). The factorisation in terms 
of the invariant /ii — 3fi2 is typical of spin ensemble averages, both in linear and mildly non-linear regime, as far as one 
investigates the rotational properties of objects with preselected inertia tensor. We remark that the lowest-order perturbative 
correction to (L'^'^) grows as D(r) E{t) oc oc t^^^ which contrasts with the Gaussian case where the growth rate is t^'^^^ 
instead (cfr. Paper I). 

Given a particular probability distribution for the linear gravitational potential tp^^\ one can compute the explicit value 
of S(i?), which involves an integral over the bispectrum. Alternatively, one can use Poisson's equation to obtain a relation 
between the potential and density bispectrum and compute the integral over the density bispectrum B^^' instead. Indeed, 
Poisson's equation implies the following relation between the two spectra 

i3l''(Pl,P2)=P?P2|pi+P2l'B^''(pi,P2) . (23) 

The latter expression is useful whenever a non-Gaussian distribution is assumed directly for the density field, rather than for 
the underlying gravitational potential field. For example, an interesting such case is the lognormal density distribution which 
was recently proposed by Coles and Jones (1991) as a reliable empirical fit for the present day density distribution function. 

Before performing the detailed calculations of the non-linear corrections to the spin, we need to describe the properties 
of those non-Gaussian fields which we intend to investigate. Unfortunately, this next section turns out to be rather long and 
readers not interested in the murky details can skip directly to subsection 3.2, where the final expressions for the spin averages 
are summarised. 
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3.1 Non-Gaussian initial conditions 

We will consider non-Gaussian statistics that are obtained from non-linear local transformations of an underlying Gaussian 
random field. These statistics are assumed either for the potential or directly for the density field and we will start with the 
first case. We restrict ourselves to the same non-Gaussian statistics numerically explored by Moscardini et al. (1991; see also 
Matarrese et al. 1991; Messina et al. 1992; Moscardini et al. 1993), namely the lognormal (hereafter LN) and the Chi-squared 
(hereafter x^) model. These statistics are imposed on the primordial potential, i.e. before the linear transfer function has 
transformed the initial scale free power spectrum to its present day shape. As an example, we will specialise to the Cold Dark 
Matter (CDM) transfer function, although our method is not restricted to CDM. 

This primordial potential is a gauge-invariant (up to a zero-mode) variable whose statistics are in principle fixed by 
primordial processes, e.g. in the inflationary scenario. This potential is not required to have zero mean, since only its gradient 
is physically meaningful. This allows for a wider variety of underlying statistics than if the probability distribution were to 
be imposed directly on the fluctuating density fleld, since the latter is required to have zero mean and in addition needs to 
be positive, p > 0, in each point. 

The models we consider are obtained by the following procedure: ip^^'' is found by convolving a stationary random fleld 
(^(q) (the primordial potential) with a real function T(q) (the transfer function): 



V'^^'(q) = J dq'T(q'-q)<^(q'). (24) 

The primordial potential ip is calculated by applying some non-linear operation on an underlying Gaussian random process 
G, e.g. if Qc for and if oc exp(G) for LN statistics. The transfer function T(q) is characterised by its Fourier transform: 

T(p) = j rfqT(q) e-^P-'' = T{p) FM''^' , (25) 

where T{p) is the CDM transfer function (e.g. Efstathiou 1990) 

T{p) = [1 + {ap + (bpf/^ + (cpf) , (26) 

with a = 6.4 (n/i^)"^ Mpc, b = 3.0 (Qh^)'^ Mpc, c = 1.7 {Q.h'^)'^ Mpc, and u = 1.13. The function F^(p) is a positive 
correction factor applied to obtain the required CDM power spectrum for t/j^^^ when the power spectrum of ip is assumed to 
be scale free: 

P^ip) = [Tip)YPM = Ap-' [T{p)Y , (27) 

and Pi>{p) ~ p~' for p 0. The latter relation also defines the normalisation constant A, which has to be appropriately 
calculated in each model since it depends on the underlying spectrum Pip{p). We will show that it is possible to factorise 
the dependence as P^{p) ~ A^F^{p) p'"^ , from which we obtain A ~ A^ F^{p)~^ F^{p) — A^, where the scale-dependence 
introduced by the function F^p cancels out. Consequently, the normalisation of i/; follows directly from the amplitude of ip. 

The assumption implicit in equation ( p^ ) is that the statistics of V''^' are of primordial origin, i.e., they refer to the 
time when the fiuctuations were outside the Hubble radius. In addition, the overall sign has to be fixed by the primordial 
physical mechanism inducing the fiuctuation field. The latter statistic is left invariant in the linear regime, independent of the 
subsequent processes which modify the wave content of ip^^^ on all scales smaller than the Hubble radius. 

Note that both and LN statistics belong to the same general class ip{ci) oc |1 -I- a~^G(q)|"; in fact, for a = 2 and 
^Q2y/2 = tJie -^2 (;jigtribution is recovered, while a ^ oo yields the lognormal one. For small wave numbers, where 

T{p) ~ 1, the gravitational potential of scale-invariant distributions obeys the scaling law: 

(i^^'Hfipi) ■ ■ ■ V;(''(mP]v)) d(MPi) ■ ■ • d(MP]v) ~ {^^'Hpi) ■ ■ ■ i'^^HPN)) dpi ■ ■ • dpN, (28) 

for any A'^, up to logarithmic corrections which arise as a consequence of the necessity to perform these integrations over a finite 
p volume. This gives a straightforward generalisation of the Harrison-Peebles-Zel'dovich scale-invariance to non-Gaussian 
statistics. 

To compute the non- linear evolution of the tidal angular momentum we need to compute the primordial bispectrum B^^ , 
which in turn is determined by the bispectrum B^, defined as usual. Since ip''^\p) — T{p) (p(p), we find a relation between 
the "transferred" and the primordial potential bispectra: 

si''(Pi,P2)=T(pi)T(p2)T(|pi+p2|)B^(pi,P2) . (29) 

It should be stressed again that the density field and the gravitational potential do not, in general, obey the same statistics. In 
addition, one cannot self-consistently assume any non-Gaussian distribution for the density field 5\ : for example, a Chi-squared 
density distribution has a positive definite 2-point correlation function, at variance with the assumed large-scale homogeneity. 
In contrast, the same statistics for ip do not lead to an inconsistency. Let us now describe the various models in more detail. 
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3.1.1 Chi-squared model 

The distribution for is obtained from the foUowing transformation of the underlying Gaussian field G: 

(^(q) = G(q)2 , (30) 

where G has zero-mean and variance of. The power spectrum of the field G is chosen as 

Pg(p) =P"'n(p) . (31) 

Here we introduced the operator n(p) = n(|p|) = i9(p — pm) ^{pm — p), where i}{x) denotes the Heaviside step function. The 
result of projecting with n(p) in equation ( pl| ) is that Pg{p) = in the interval \pm, Pm] and zero outside. The cutoff 
constants Pm and Pm are introduced to avoid divergences at both infra-red and ultra-violet 'wavelengths' . In terms of these 
cutoffs, the field G has variance cr^ = {2-K'^)~^\n{pM /Pm). On the other hand, the constant (fio has the same dimensions 
as if and can be either negative (x^ model) or positive (xp model). These names are chosen according to the sign of the 
corresponding skewness (5^) /{S'^)^^^: the Xn iXp) model is negatively (positively) skewed. Finally, it is worth noting that the 
fiicker-noise choice for the Pa spectrum, giving equipartition in Fourier space, is stable up to negligible corrections against 
the non- linear transformation originating ip, both in Chi-square and LN non-Gaussian models (see below). 

The Chi-squared model is an example of a scale-invariant non-Gaussian statistic (Otto et al. 1986; Lucchin & Matarrese 
1988). Such models have been considered by Coles & Barrow (1987) in connection with the two-dimensional distribution of 
GBR fluctuations. In the inflationary context, Bardeen (1980) proposed a model where adiabatic perturbations are described 
by a squared Gaussian process, our Xp, whose G fleld has a non-scale-free power spectrum. 

It can be shown that the ip power spectrum corresponding to the choice ( ^l| ) is (see Appendix B for details) 

PAp)^^ [/3(p) + 21n(l+p/p„)-l]p-«n(p), (32) 

where /3(p) = (l — p/pm)^ , for pm < p < 2pm, and (3{p) — 1 + 2 ln(— 1 +p/pm), for 2p„i < p <C Pm- The normalisation constant 
A is related to ipo as j4 « — <^o/27r^ . 

Relevant to our investigation is the expression for the linear gravitational bispectrum B^, which turns out to be given by 
a convolution over the underlying spectrum Pq, namely (see Appendix C for details) 

5^(Pl,P2) = || / dpPG{p)PG{\pi+p\)PG{\p2-p\) . (33) 

Note that the statistics of the Chi-squared potential ip are independent of the variance of the underlying Gaussian field 
(see Appendix C). This is not the case for the LN distribution (see below), for which the departure from Gaussianity can be 
tuned by varying . 

Finally, we define for later use the operator 

n(pi, P2) = n(pi) n(p2) n(pi + pa) , (34) 

i.e., the projection n(pi, P2) is zero if any of the two vectors pi or p2 or their sum has modulus < p™ or > pM and is one 
otherwise. 



3.1.2 Lognormal model 

The LN distribution for ip> is obtained from the following transformation of the underlying Gaussian field G: 

^(q) = <^oe«('" , (35) 

where, as before, G has zero-mean and variance a^. Note that (ip) = fo exp(cif /2) and there is again the freedom of assuming 
ipo to be negative (LN„ model) or positive (LNp model). 

The lognormal statistics has been analysed in the cosmological context by many authors (see e.g.. Coles 1989; Coles & 
Jones 1990; Messina et al. 1990; Coles , Melott & Shandarin 1993; Catelan et al. 1994; Sheth 1995). 

In contrast to the case, one gets lognormal statistics which differ by varying amounts from the Gaussian distribution 
by tuning the variance of G: for small a^, the non-Gaussian character of ip> is manifest only in the properties of rare high 
peaks, whereas for larger cf,, the power spectrum P^ deviates strongly from fiicker-noise. 

Next, we analyse the hierarchical structure of the correlation functions of the field ip. The 2-point correlation function 
{V3(qi) V3(q2)) - (vs)^ = iAl) is given by 

Uq) = (^)' [e^"^'"' - 1] , (36) 

where ^0(9) = (G(qi) G(q2)) and g = |qi — q2|. The 3-point correlation function {{ifi — {if)) {(p2 — {if)) {ipz — {if))) = C<^(123) 
is (see Appendix C for the derivation) 
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C^(123) = [U^2) ^^(13) + e^(12) C^(23) + ^^(13) C^(23)] + {cp)-'' [U^2) C^(13) C^(23)] 

The power spectrum Pip{p) is given by 



1] 



(37) 



(38) 



and the formal expression for the bispectrum is 

{<^)='B^(pi, P2) = {^f {PM PM + PA\pi + P2I) [pApi) 



pAp2)]}+ J 



dp 



{27V) 



^ PAp) PfilPi + P\) PA\P2 ' P\) ■ (39) 



A remark is now appropriate. Since we are interested in comparing the results of these non-Gaussian investigation against the 
results from Gaussian initial conditions reported in Paper I, it is important to preserve the form of the power spectrum. On 
the other hand, it can be easily understood from equation ( ^s]) that the spectrum of the field is a complicated function of 
the wavevector p, even for the simple underlying spectrum Pg{p) adopted in equation (|3l|). Fortunately, the transformed field 
ip may be easily constrained to have the required power spectrum if we restrict ourselves to those LN distributions for which 
of < 1. This restriction leads to some loss of generality, but the formalism becomes much simpler, and, above all, the power 
spectrum of the underlying Gaussian field is preserved. In this case, since iail) is a decreasing function of q, the condition 
cr? < 1 implies that |^G(g)| < 1, and a Taylor expansion of the exponential in the integrand ( ^ ) leads to P^(p) ~ {ifi}^ Pg{p) ■ 

1. Consequently, the power spectrum may be written 



In the following we restrict ourselves to LN distributions with 
as 

pAp) ~ w 7 ^ p' n P . 

and the bispectrum as 



2 < 



(40) 



-B<^(Pl,P2) 



27VW 



ln{pM/Pm) 



dp 

ln(pM/Pm) J {2ny 



-3 -3 , 
Pi P2 + iPi 



P2I 



-3/ -3 I -3\ 
[Pl +P2 ) 



\pi+p\ 1p2-p1 'n(pi,p)n(p2,-p) 



n(pi,P2) 



(41) 



Note that the bispectrum _B^(pi, P2) is zero outside the region pm < Pi,2 < Pm and pm < |pi + P2I < Pm- The potential ip 
now has the CDM power spectrum with amplitude A ~ 2-K^a^ (pi e°^G / \n{pM /Pm)- 



3.1.3 "Linearised" lognormal distribution for the density field 

A completely different approach may be adopted if the non-Gaussian distribution is imposed directly on the linear density 
fluctuation field 5\, rather than on the underlying potential. Although one is not allowed to choose any non-Gaussian density 
distribution, as argued earlier, there are some advantages in considering models where non-Gaussianity is coded directly in 
the statistics of the matter distribution. An interesting case is the LN density distribution, proposed by Coles & Jones (1991) 
as a reliable empirical fit to the present day density distribution function. A theoretical explanation for this fact followed 
from the analysis of Bernardeau (1994), who computed the evolution of the one-point probability distribution function of the 
cosmic smoothed density induced by gravitational evolution. Indeed, his results indicate that, during the mildly non-linear 
regime, the density distribution is fairly well approximated by a LN statistic when the spectral index n is close to —1. For 
larger values of n or for smaller values of the rms cr^, the resulting LN distribution is close to Gaussian. The formation of 
large-scale structure by gravitational instability from primordial lognormal density fields has been analysed by Weinberg and 
Cole (1992) using A^-body simulations. A further reason to consider a LN density field is motivated by the recent extension of 
the Hoffman-Ribak algorithm (1991) for the construction of constrained Gaussian random fields by Sheth (1995) to include 
also lognormal statistics. 

Let us analyse more closely the complex hierarchical structure of the LN density distribution. In particular we wish 
to discuss the approximate expression for the 3-point LN correlation function we are allowed to adopt in our investigation. 
Fortunately, as we will see below, the lognormal distribution permits the constraining of the deviations of p[5) from the 
Gaussian distribution uniquely in terms of the linear-regime condition. Let us suppose that the (linear) density field p = p(,[l+(5] 
is obtained by the transformation 

p = P6e^-"G/\ (42) 

where, as before G is an adimensional zero-mean Gaussian random field with variance a^. The mean value is given by (p) = ph. 
We stress the fact that equation ( ^ ) is mathematically identical to the transformation (^jj), but, since it involves the density 
field, it is physically very different. 
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The 2-point LN correlation function ^^(g) is given by 
6(g) = (5i(qi)<5i(q2))=e««(')-l, (43) 
where ^g(O) = a^. Evaluating ^s{q) at zero lag gives the variance of the LN random field 61: 

aj = e'b' - 1 . (44) 

We stress that the equations (^^ and (^) do not contain any time dependence and combine two concepts, that of "non-linear 
regime" and that of "deviation from Gaussianity" . Since for a LN distribution the A^'-point (TV > 2) functions can be expressed 
as a Kirkwood expansion of the 2-point functions (see Coles & Jones 1991; also below), the linear-regime condition can be cast 
simply in the form (t| ^ 1 which then also implies <C 1, i.e. the primordial LN density distribution with small fluctuations 
cannot deviate strongly from Gaussianity. Alternatively, if the primordial LN density distribution were to deviate strongly 
from the Gaussian one then the initial conditions would contain an intrinsic non-linear signal, unacceptable on observational 
grounds. Finally, since £,G{q) is a decreasing function of g, c[? ^ 1 implies |^G(g)| <^ 1 and consequently also |^i(g)| <C 1, 
which is important for the subsequent discussion. 

We first show that the density power spectrum and the underlying Gaussian spectrum Pg(p) are similar in the approxi- 
mation (Jq ^ 1. Indeed: 

P6{p)~Pg{p) = Ap[T{p)]'iI{p), (45) 

where T{p) is the CDM transfer function. Additionally, we derive some properties of the structure of the LN linear potential 
bispectrum B^, which is the Fourier transform of the 3-point correlation function C« (lii I2, qa) = ^5(123), 

Ca(123) = is{l2) C5(13) + ^^(12) ^^(23) + £,s 5^(23) + ^a(12) ^^(13) C« (23) . (46) 

The important point to note is that the expressions (^sj) and ( ^ ) may be simplified if we restrict ourselves to the case cr^ ^ 1, 
which corresponds to the linear regime, as explained earlier. Since in this case also <^ 1, it seems reasonable to approximate 
the 3-point function C,s by the relation 

a(123) ^^^(12)5^(13) +6(12)^^(23) +5^(13) C^(23) . (47) 
The corresponding bispectrum reduces to 

B,(pi,P2) ~ Pi(pi)P5(p2)+P^(|pi +P2I) [P^(P1)+P^(P2)] 

« Pg{pi) Pg(P2) + Pg{\Pi + P2I) [Pg{pi) + Pg(P2)] , (48) 

where the cubic terms have been neglected with respect to the quadratic ones. This approximation allows a considerable 
simplification of the calculations, yet the loss of accuracy of the whole description is small. We stress that the same strategy 
cannot be applied to the LN gravitational potential since the constraint a^p ^ 1 does not imply a restriction to the linear 
regime. 



3.2 Ensemble averages 

3.2.1 Linear approximation: (L'^-*^) 

The ensemble average of the rms linear angular momentum, (L'^'^), is computed and discussed extensively in Catelan & 
Theuns (1996a). Here, we briefiy summarise the main results, which hold for any underlying probability distribution, i.e. for 
both Gaussian and non-Gaussian statistics. We have: 

(L(i'2) = Ai)(r)2(M?-3M2)a(7?)\ (49) 

where the quantity <j{R)^ is the mass variance on the scale R which is obtained from the power spectrum through o-{R)^ = 
(27r^)"^ dpp^ P^{p)W{pRf . We wiU adopt a Gaussian smoothing function, W{pR) = exp(-p^_R^/2). Equation (^) holds 
for any power spectrum {^'""^'(pi) ip^^'' iP2)} 4, = (27r)^ (5d(pi+P2) P^(pi) and the value of (L'"*^'^) depends on the normalisation 
of the spectrum. The general expression (^) is independent of the details of the boundary surface of the volume F but depends 
on fii and fi2, which are the first and the second invariant of the inertia tensor Jap. Specifically, denoting the eigenvalues of 
the inertia tensor by ti, 12 and ta, one has /ii = ti + 12 + ts and fi2 = iiL2 + tits + 12^3 . For a spherical volume, li — L2 ~ I'S, 
hence (L'^-''^) cc pl\ — 3/i2 = 0, as we stressed before. 



3.2.2 Higher-order approximation: (L'""^' • L'^^) 

The calculation of the term (L'^'.L'^') takes advantage of the results of the second-order approximation, and essentially reduces 
to specialising the function E(_R) to the particular statistics chosen, since E(i?) = IhlC^xyxy ffi as explained in the intro- 

duction of Section 3. The final expressions, valid for a CDM power spectrum Ps{p) = Ap[T{p)]'^ — Ptpip) ~ P* [^(p)]^ PfiP)^ 
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where T(p) = F^(p) ^^'^ T(p) and P<fi(p) ~ F;p{p)p , may be written for the different statistics as (compare with equa- 
tion!^): 

• For a Chi-squared gravitational field: 

(L(i) . L(^') = A D{r)Eir) - 8^2) E^. (R) , (50) 
where the function E^2 depends on the smoothing scale R and the normalisation of the spectrum A, as 
E^2(i?) = 15/Ci'j,>,j,(pi,p2)eT(pi)T(p2)T(|pi +p2|)B^(pi,P2) . (51) 
The correction F^{p) is such that [see equation ( ^^ for the symbols] 

T{p) = T{p) [P{p) + 2 ln(l + p/p„0 - 1] , (52) 
and the normalisation constant is A = </?o/2vr^. The primordial bispectrum may be deduced from equation (p3[), 

s^(pi,P2) = ±(2yi)'/' / dpp-'|pi + pr^^|p2-pr'n(pi,p)n(p2,-p)n(pi,p2) . (53) 



Typically, these integrals have to be evaluated numerically. The upper (lower) sign is for positively (negatively) skewed dis- 
tributions. Note that the average (L*-^' • L'^'^-') factorises the same invariant /if — 3fi2 of the inertia tensor J as appeared in 
the linear term [see equation (^^]: this is a direct consequence of the isotropy of the universe which is independent of the 
underlying statistical properties of the matter distribution. 



• For a LN gravitational potential field, the calculation may be partially done analytically to give E_ljv(-R), which is a sum of 
three terms: T,ln (R) = E^"j^(i?) + 'S'^llfiR) + E*J'],(i?) which originate from the Kirkwood expansion of the bispectrum B\ 



m equation (Ba): 

Eijv(ii) = 15K:^^ly{pi,p2) e r(pi) T(p2) T(|pi + pal) B^(pi, P2) ; (54) 
where 



B,(p„P2) - ±A^^--I ^^''""o 



In (pAf/Pm) 



P1'P2' + |P1+P2| '(Pl'+P2') 



+ T"" ^ li^P 'Ipi+Pl '|P2-pi 'n(pi,p)n(p2,-p) 

ln(pM/pm)j {2tiY 



n(pi,P2). (55) 



In this expression, the product p^ '^pj ^ gives rise to the E^°)^(_R) term, E^''^(J?) corresponds to ~ |pi -I- pal ^{p^ ^ + P2 
and finally E^'^^(i?) is the loop contribution. In Appendix D we will show how these integrals may be partially done analytically. 

• For a "linearised" LN density field, we obtain in a similar way as before the function T,s{R) ~ E^"-'(ii) -I- E^'(i?), which is 
the sum of two contributions which originate from the two terms in the linear bispectrum Bg^^ in equation (^^ : 

Ss{R) = 15/Ciy'^y(pi,p2) ep^Vj^lPi +P2r^-Bf'(pi,P2) , (56) 
where 



^1" (Pl, P2) = A^ n(pi, P2) {pi P2 T(pif T{p2f + |pi + P2 I T(|pi + P2 I 



Pir(pi)'+P2r(P2)' 



(57) 



Again, E^"' corresponds to the pip2T{pi)^T{p2)^ term and E^**' is the remaining contribution. Also these integrals may be 
reduced to simpler expressions analytically: the final formulae can be found in Appendix D. 

Numerical values for the various Es are shown in Fig. |l|. Only positively skewed contributions are shown, negatively skewed 
have opposite sign. Some terms are logarithmically divergent for pm —> 0, i.e., for large scales. This divergence can be traced 
back to the bispectrum. Note that the limits pm < Pi,2 < Pm follow from restricting the initial scale-free power spectrum Pq 
to a finite interval; however, this generally sets no restrictions on |pi -I- pal and the divergence occurs for elongated triangles 
Pl —p2 which can be controlled by the infra-red cut-off Pm in n(pi, pa). There are no ultra-violet divergences {pm — > 00) 
because such small scale interactions have been explicitly smoothed over. The different Es have been scaled by appropriate 
powers of the mass variance o-{R) to make them independent of the amplitude A of the power spectrum. In addition, the 
LN corrections have been scaled by their dependence on the dispersion c;-, of the underlying Gaussian model. Finally, some 
Es have been scaled by powers of R to reduce the scale dependence of the plotted quantity. The E^'^]^ and E^2 contribution 
follow from 6D Monte-Carlo integrations and are less accurate than the others, which involve 3D numerical integrations. 
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Figure 1. Correction terms to the linear angular momentum versus scale R, for several non-Gaussian statistics indicated in the panels 
(only the positively skewed ones are shown). The E-terms have been scaled by powers of the mass- variance cr(R) to make them independent 
of the spectral normalisation and by powers of R to reduce the scale dependence of the plotted quantity. The lognormal terms have in 
addition been normalised by powers of c^. Full lines corresponds to the choice pm = 0.005Mpc~^ and dashed lines to pm = O.OlMpc"^; 
Pm = 2.5Mpc~-'^ for both, the scale R is in Mpc {h = 0.5). 
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Figure 2. Left panel, lower scale: skewness 5 as a function of scale R (in Mpc, h = 0.5) for the different (positively skewed) statistics: 
Sln and Ss denote the density skewness for the lognormal gravitational potential and the lognormal density field, respectively. The 
individual contributions to S^jf are scaled by powers of q~., as indicated; Ss is scaled by a. The dotted line (upper scale) is Ss versus 
cr/cg for the LN density field taken from equation (43a) in Bernardeau & Kofman (1995). Filled diamonds show our determination of 
S^"' + S^''^ versus cr/ag for comparison. Right panels, lower scale: skewness S (upper panel) and kurtosis K (lower panel) for the 
(short dashed) and LN (long dashed) unfiltered models as function of . The values are independent of . The filtered skewness for 
the LN gravitational potential distribution is plotted for comparison (upper panel, full lines) for filtering scales (from top to bottom) 1, 5, 
10 and 50 Mpc, with h = 0.5. The dotted lines (upper scale) show the functions in equations (43a) and (43b) from Bernardeau & Kofman 
(1995) for skewness (upper panel) and kurtosis (lower panel) versus cr/ag for the LN density field. For all curves, pm = 0.05Mpc~'^, 
Pm = 2.5Mpc-i, h = 0.5. 
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4 ANGULAR MOMENTUM AT THE MAXIMUM EXPANSION TIME 

In this section we quantify the relative non-Gaussian perturbative corrections to the linear angular momentum for the different 
statistics by computing the parameter 

at the maximum expansion time, D(tii)g[R) — 1 (we computed the maximum expansion time for the spherical model; see 
Peebles 1969; see also Paper I). 



• For the Chi-squared distribution we obtain 

T,.(i?) = II [Dir^MR)] ^ ~ ±2.1 (R/h-'Mpc)-'^ , (59) 



where the last value is estimated from Fig. (|l|). This correction is of 0(1) on galactic scales. We suggest that this indicates 
that higher-order terms in the perturbative expansion probably contribute significantly. We reiterate that the upper (lower) 
sign refers to positively (negatively) skewed distributions. 

• For the LN distribution we obtain 

Ti^(i?) = ^ ^itiR) + ^tUR) + ^^lUR) [D{rn)a{R)] ~ ±[OAa^ + a^{R/h-'Mpcy/^ + 0.2a^{R/h-'Mpcr% (60) 
DD ^\R) 

which depends explicitly on q-,, i.e., on the level of non-Gaussianity imposed on the initial state. Below we estimate an upper 
limit for ac and hence an upper limit to this correction. 

• For the linearised LN density field, the scaling is rather different. This is because the degree of non-Gaussianity is determined 
directly by the amplitude of the power-spectrum and consequently the non-Gaussian correction term must also depend on 
that amplitude, parametrised for example by the mass-variance in spheres of Rg = 8/i~^Mpc at the present time, erg. This 
complicates the dependence of the correction on the mass scale M. Approximating the scaling of a with R for the standard 
CDM model as a{R)/aa « (_R8/_R)i i+0 25i°g(«/«s) ^^d given the fact that the relevant Es scale cx Ra{R)'^, we obtain 

2E E^;\R) + Ef\R) 2E sW(fl) + st'')(i?) . y. 1+0.25 logWi^e) 
= [DimMR)]^^ ^-^^ Ras ^-j ~±6.3as, (61) 

where the numerical value is appropriate for scales R = 0.5^~^ Mpc. Clearly, LN non-Gaussianity imposed directly on the 
density field can change quite drastically the amount of angular momentum induced on galactic scales. The absolute value of 
the correction depends on the normalisation of the power spectrum, through as , which contrasts with the other distributions 
which are independent of erg. 



4.1 Estimating the degree of non-Gaussianity 

The relative spin correction Tln for the case of a lognormally distributed gravitational potential ip depends on the free 
parameter ac, cfr equation (|6o|), which fixes the degree of non-Gaussianity. An upper limit on the allowed degree of non- 
Gaussianity, and hence on ac, can be obtained by analysing the deviations from Gaussianity for the density field. This also 
allows us to show that distributions with similar skewness and kurtosis can nevertheless induce widely different corrections 
to the angular momentum. 

In this spirit, we have computed numerically the skewness S = {5^)/(5^)^^^ as a function of scale for all the considered 
non-Gaussian models starting from their respective bispectra: the result is shown in Fig. (^) which also shows the unfiltered 
kurtosis K = ((5*)/(5^)^ — 3. The dotted lines denote the expressions given by Bernardeau & Kofman (1995) for unfiltered 
skewness and kurtosis versus a for the LN density field (their equations (43a) and (43b)). From Fig. (left panel) it is clear 
that the skewness decreases weakly with scale for all models, except for the "linearised" LN density field which decreases 
more strongly: for a ^ 1, Sln ~ 3cr (and Kln ~ 16(j^). Note that our estimate for the skewness for the filtered "linearised" 
LN density model versus a agrees well with the value for the unfiltered one. 

We can estimate an upper limit to for the gravitational potential LN model by requiring that the primordial skewness 
of the density field, 5*,^^^ , is much smaller than the observed skewness, S ~ 4 on scales J? ~ 8 Mpc. We thus set 

S^^Jq,)^Aa, (62) 

with a <C 1. Using the approximations Sy,^^ [a^] ~ l.Sq, -I- O.Ga^ (from Fig. 2, for ^ 1) we get 

cTg « 2.22a - 3.65a^ + 18a^ for a«l . (63) 
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Taking a ~ 0.1, we obtain ~ 0.2 from which TiAf(_R = 0.5h ^Mpc) ~ ±0.33. We conclude that in this approximation the 
non-Gaussian correction to the linear angular momentum for a LN potential is of the order of < 33 per cent. 

Summarising, it is clear that different underlying statistics give rise to differing corrections to the linear estimate of 
the angular momentum. For instance, the contribution of higher-order spin corrections appears to be non-negligible for 
and for the "linearised" LN density field, whereas the porturbative expansion for the LN potential field scorns to converge. 
Clearly, more detailed investigations of these non-Gaussian corrections seem warranted. Finally, note that in all cases positively 
skewed distributions increase the angular momentum of collapsing objects with respect to the linear estimate, and vice versa 
for negatively skewed ones. 

5 SUMMARY AND CONCLUSIONS 

In this paper we analysed the corrections to the linear growth of the tidal angular momentum L acquired by a proto-object 
(protogalaxy or protocluster) as a consequence of non-Gaussian initial conditions. We have used Lagrangian perturbation 
theory for the displacement field to obtain perturbative corrections to the linear angular momentum. Whereas the linear rms 
angular momentum involves integrals over the power spectrum alone and hence is independent of the assumed statistics of 
the random field - the lowest-order perturbative corrections involve integrals over the bispectrum and so do depend on the 
underlying probability distribution. We showed that for a generic non-Gaussian distribution the corrections to the variance 
of the angular momentum grows oc t^^^ , which contrasts with the Gaussian case for which the growth rate is t^'^^^ (both for 
the Einstein-de Sitter universe). This is a consequence of the fact that the lowest order perturbative spin contribution in the 
non-Gaussian case arises from the third moment of the underlying density field, which is identically zero for a Gaussian field. 
In this formalism, the resulting expression for the variance of the spin factoriscs the same shape parameter as in the Gaussian 
case, smce oc /jf — 3/J2, where /Ji and are the first and second invariants of the inertia tensor of the collapsing 

object. This result does not depend on the assumed statistics, since it is a consequence of the isotropy of the Universe (see 
Appendix A) . 

The spin variance was evaluated explicitly £is a function of scale for a variety of multiplicative non-Gaussian statistics, 
namely and lognormal distributions imposed on the gravitational potential and "linearised" lognormal statistics assumed 
directly for the density field. We characterised these corrections in terms of T(_R) = (L^)/{L''^-' ) — 1, i.e., in terms of the 
relative contribution of non-Gaussian corrections to the linear spin. In general, during the mildly non-linear regime, positively 
skewed distributions increase the angular momentum with respect to the linear term, and vice versa for negatively skewed ones. 
However, the convergence properties of the perturbative series depends strongly on the assumed statistics. We found that the 
and "linearised" lognormal density statistics induce 0(1) corrections to the linear angular momentum, presumably indicating 
that higher-order terms, not discussed in this paper, are non-negligible. Consequently, angular momentum acquisition in 
models with this kind of underlying statistic appear analytically intractable. In contrast, the investigation of the lognormal 
distributions imposed on the gravitational potential seems more promising in this respect, since the lowest-order non-Gaussian 
spin correction are of order 30 per cent, suggesting convergence of the perturbative series. The calculations in this paper do 
not take into account non-gravitational processes or highly non- linear events which may play important roles in the formation 
of galaxies. 
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APPENDIX A 

In this first appendix we report the details of the derivation of equation for the spin correction (L*^^-*- L'^'), namely 
(L(^) . L(2)) = A D{t)E{t) (m? - 3^2) , (64) 
where the wave-vector contribution E(_R) is defined by 

E(i?) = 15 (Oi'j,) V'i^) = 15 K.'ily e B*,'^ . (65) 
We will show that this expression is a direct consequence of the isotropy of the universe. 



Our starting point is the general expression [see equation (|18[)] 
(L^^) . L'^') = b{r) E{r) e^p, e^p,,, J^, J„,y {0<t' T^'-^l) ■ (66) 
The Levi-Civita tensors may be written as 6^/37 ea/3'7' = SpfjiS^-yi — Sp^iSfji^, which leads to: 

^L(i) . L{2)^ ^ ^(^) Eir) J^-, [J^,, {Vf^ Vf^,} - J^,, {Vf^ O^^l)] . (67) 

Without loss of generality, we now assume that the chosen frame of reference coincides with the eigenframe of the inertia 
tensor, in which case Jafs — LaSap, where La are the eigenvalues of J (no summation over a in this specific case!). In this 
particular reference frame, we get 

^L(i).L{2)) ^ I)(r)i5(r)(.^-...,){0(t'2'^l'> (/? / 

= D{t) E(t) [2(4 + 4 + 4) {Vg^ ©(2)) - 2{l,^Ly + L^L, + LyL,) 1,(2))] 

= 2DiT)E{T){pl-3fi2){Vi'-Jvi^J) , (68) 

where ni = Lx + Ly + and /12 = tajty + l^Lz + Lyi^ are the first and the second invariant of the inertia tensor J . The latter 
equation (168) was used in Section 3.2 



APPENDIX B 

In this Appendix we describe how to obtain the expression of the Chi-squared power spectrum in equation (B2), using the 
same notation as in Section 3.1. The power spectrum of the 2-point correlation function of the field ip = ipo G^s given by 

„2 



P^{p) = 2vi I dq[^Giq)r e-^''-'' ^ ^ / dk Pa (k) Pg {\k - p\) 



(69) 



Changing variable from azimuthal angle (using spherical coordinates in k) to £ = |k — p| simplifies the previous integral to: 



dkkPaik) 



p+fc 



dllPoil) 



(70) 



Let us now consider the fiicker-noise power spectrum for G defined in equation (^ij) and introduce the parameters e = km/p 
and Q = kulv- The condition kra < p < ku then translates into e < 1 and Q > 1 and the previous integral simplifies further 
to 



27r2p3 



dk k 



Min{Q,fc+l} 



Max{e,|fc-l|} 



In the limit Q ^ 00 one has 



PAp) 



27r2p3 



dk 



+ 



P Max{e, 1 - fc} 7^ fc2 Max{e, fc - 1} 



dk 



dk 1 

p FT 



and two expressions for P^ may be recovered: for p > 2fcm, i.e. e < 1/2, 



P^P) = 



2n^ p3 



2 In 



while for km < p < 2fcm, i.e. e > 1/2, the power spectrum reduces to 



pAp) = 



27r2p3 



2 In 



+ 



1 - e 



(71) 



(72) 



(73) 



(74) 



Substituting e = km/p into the latter equations leads to the compact expression for the Chi-squared power spectrum in 
equation (|32|). 
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APPENDIX C 

In this third Appendix we summarise those properties of the Chi-squared and lognormal distributions which are relevant 
to our analysis. To indicate the non-Gaussian field, we employ the same symbol (fi used in the main text for the primordial 
gravitational potential; however, the results derived below hold for any non-Gaussian Chi-squared or lognormal field. The 
reader interested in a more comprehensive treatment is addressed to Kendall & Stuart (1977). 



C.l Chi-squared distribution 

Let G be a Gaussian variable with zero mean and variance cr?. The variable = ipoG^ has a Chi-squared distribution, with 
characteristic function 

= I ■ (75) 

The moments of tp about the origin ip = Q can be obtained by evaluating successive derivatives of the characteristic function 
$^(t) a.tt = Q: 



dt 

Explicitly, the first few moments are found to be: 



(76) 

i=0 



^ii = ifiocr^ ; ^i2=3ifla^; fj,3 = ^5 ipl a^' . (77) 
The reduced moments about the mean, or cumulants, are {ki = 0) 

K2 = {{ip - {(p)f) = fi2 - = 2 ifl , (78) 

'^a = {{f - i'fi})'^} = - S fJ'i tJ-2 + 2 = 8 (fil , (79) 

The cumulants k„ with n > 2 vanish, by definition, for a Gaussian distribution. The deviation from non-Gaussianity of a field 
ip with given dispersion cr^ obtained in this way cannot be tuned by changing of the underlying Gaussian field. Indeed, 
since ct^ — Stp^a^ is assumed to be constant, we find K2 = 2a^/3 and K3 — 8(j^/3, independent of q,. Similar relations hold 
true for all higher moments as well. The same remark applies to all fields (x G", obtained from powers of an underlying 
Gaussian field. In contrast, we will show below that for LN statistics, the deviation from Gaussianity can be tuned by varying 

The characteristic function of if is useful to calculate the moments of its one-point PDF. To compute the correlation 
functions of the field (p, it is convenient to work with the characteristic function "1>[^' of the multivariate Gaussian probability 
distribution of the field G(q), where the dependence of G on the 'spatial' variable q is now emphasized. It is given by: 



<I.(,"'(t)=exp 



2 



(80) 



where t = (ti, . . . ,t„), t denotes the transpose of t and C is the covariance matrix of the Gaussian field G, whose elements 
are 

{G(qOG(q,)) =eG(g,,) (81) 

and qij = | qi — qj | . Here, (g) is the covariance (correlation) function of the Gaussian field G with (0) — . By definition, 
the 2-point correlation function of the field ip is 

U<l) = {{fi - if)) i'P2 - (p))) = {pi P2) - {pf , (82) 
where pi = p{cii)- We have 

^-^{^1^2) = {G?G^)= 1^ ^a^ + 2^G{qr. (83) 
Combining equations (p2[) and (p3|), we obtain the following relation between the correlation function of p and that of G: 



Wl) = 2p,i^G{qr ■ (84) 

In a similar way, using the characteristic function of order n of the multivariate Gaussian distribution function, one may 
derive the higher-order correlation functions of the field p. Thus, the 3-point connected correlation function is 

C^(123) = {{pi - {p}) {p2 - {p}) ips - {p))) = 8^Ug(12)Cg(13)Cg(23) , (85) 



where Cg(*j) = ^G{<ii,<lj), and so on. Note that, from equations ( |84| ) and (|85[), one has K2 — £,ip{0) ~ 2p^a^ and 
K3 — C'p(O) = 8po cr?, in agreement with equations ( [t^ and (^^. The bispectrum may be recovered by performing 
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the Fourier transform of equation (pq), which gives the expression reported in equation (p^). We stress once more the absence 
of any exphcit dependence on the underlying variance a^. 

C.2 Lognormal distribution 

Let us again suppose that G is a normal variable with zero mean and variance rr^. The variable ip = ifioe'^ is said to be 
lognormal distributed. The 2-point correlation function of the field ip, defined in equation (p2|) , may be obtained from the 
correlation of the field G via the bivariate characteristic From equations (^) and we get 

(^-2{VPiV92) = {e^i+^^) = -I>^''(-i,-i)=e€e««('') . (86) 

Moreover, (p^^ {<p) = i) = e'h^^. Combining these expressions one gets again a relation between the two 2-point 

correlation functions: 

U<1) = ivf b^"^'^ - 1] • (87) 
In a similar fashion, one can obtain the 3-point correlation function, defined in equation (^^. Since 

Cv(123) = {ipi tp2 ips) - (<y3) (<y3l 'P2) - {f) {p>2 'fis) - {f) {<fl ifs) + 2 {iff , (88) 

one obtains 

{^i^2V3) = ivf + {<p) K^(12) + 5^(13) + C^(23)] + C^(123) . (89) 
On the other hand, the left hand side of the latter equation may be written as 

(.^1^2 ^3) = ^?(e«i+«^+«^) = ¥'^$(?)(_i,_i,^i) 

= (^)3e«G(12)+CG(13)+eG(23) ^ ^^^3 ^ (^^y\.^[l2)] [1 + (^)-^C^(13)] [1 + {^)-'e^(23)] 

= {^)-^ ^^(12) ^^(13) e^(23) + if,)-^ [e^(I2) C^(13) + e^(I2) C^(23) + e^(I3) e^(23)] 

+ {^) K^(12)+ 5^(13) +5^(23)] + {^)^ (90) 
and, from (^), one recovers a Kirkwood-like expansion for the 3-point connected correlation function: 

C^(123) = (if)-' K^(13) e^(23) + C^(12) ^^(13) + C,,(12) C,,(23)] + {<^)-=' [C^(12) ^^(13) ^^(23)] . (91) 
The first few cumulants k„ of this lognormal distribution are: 

AC2 = e^(0) = fU^G (e"G - l) , (92) 
«3 = C^(0) = ipl e'^'fe'/^ (e€ - l) ' (e'^G + 2) . (93) 



Taking again cr^ — K2 constant by setting ip^ = K2/e'fe (e'fe — 1) where is now a free parameter, one can write K3 = 
Kj''^ (e'c -|- l)(e'G — 1)^^'^, from which it is clear that the skewness of ip can be tuned independently of its variance. Clearly, 
this is also true for higher moments. In other words, the deviation from Gaussianity of the LN tp can be tuned by changing q,. 
In particular, in the limit cr? — ^ the lognormal distribution tends to a Gaussian. The bispectrum of the lognormal variable 
tp may be derived by Fourier transforming the expression ( ^l[ ) (Catelan, Chodorowski, Lokas and Moscardini 1997). The final 
result is reported in equation ( p9| ) . An explicit dependence on the underlying variance cr^ can be noted in equation ( pl| ) , which 
holds in the regime cr? ^ 1. 



APPENDIX D 

In this last appendix we want to show how to perform analytically some of the integrations when computing the spin correction 
E(ii) for the different non-Gaussian fields we have considered. The compact expression for E(7i) in terms of the kernel K^yxy 
enables us to reduce easily the number of numerical integrations by exploiting rotational invariance of the integrals. 

• We start by considering the simplest case, the 'Eln{R) correction for the LN gravitational potential (p: 

Eln{R) = 15 IC^ily e B*,'' = Si^ji {R) + T.%{R) + E^^ (i?) . (94) 

As discussed in the main text, the three contributions T,^^^{R) originate from the Kirkwood relation of the LN bispectrum 
B^-* [see equations ( p9| ) and (^)]. Let us consider the first of these three contributions: 

^it{R) = 15 /ci',U(pi,P2) e T(pi)T(p2) T(|pi + P2I) P^(pi) P^(P2) . (95) 
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which can be written explicitly as: 



15 /■dpidp2 (Pi +P2)^ (pi + P2; 



y ^(2) 



|Pl +P2l^ 



'^^''(Pi, P2) [W{\t>i+V2\R)Y Tipl) T(p2) T(|pi +P2I) P^(pi) P^(p2) .(96) 



We now use the relation P^{p) ~ A,^F^{p)p '^n(p), with = 1 for the LN case (and hence T{p) = T{p)) and in addition 
— 27r^cr^ {(f)^ / ln{pM /pm) ~ A. Furthermore, without loss of generality, we assume that the vector pi is chosen along the 



2-direction, in which case the previous expression simplifies to: 



(97) 



Ei^liR) = ^i5f^ [ ^^ /^f^",, ^'^'(pi,P2)[iy(|pi+P2|i?)]'r(pi)r(p2)r(|pi + p2|)pr^P2-^n(pi,p2) 

(W J l^Ti") !Pi + pal 

Note that the integrand does not depend on the azimuthal angle 9i nor on the precessional angle (pi; in addition, the <^2- 
integration may be done analytically. The final result may be cast in the form: 



15 



8(27r) 



±A 



3/2 



ln(pM/p rr 



dpiT{pi) / dp2plT{p2) / dSa (sin6l2)^ 



X (pipa' +2cos6»2) 't(|pi +P2I) [l4"(|pi +P2I-R)] n(pi,p2) , (98) 

where 62 is the angle between pi and p2. The upper (lower) sign is for positively (negatively) skewed LN distributions. The 
remaning integral needs to be computed numerically. 



The other two contributions may be calculated in a similar way. The second contribution is 
15 



4(27r) 



ln(pA//p™ 



dpiPi"'/'r(pi) / dp2p"'"r(p2) / de2(sme2y 



PiP2^ +P2Pi^ + 2cose2) ^^""rdpi +P2I) [VF(|pi + P2|i?)]''n(pi,p2) 



(99) 



The computation of the third contribution is the more involved, since it corresponds to a 9 dimensional integration. Recalling 
that the bispectrum B^p depends only on the moduli of the vectors pi and p2 and the mutual angle 9\2 ~ 62, one has 



15 



8(27r)7 



±A 



3/2 



\ ln(pA//Pm 



3/2 



dpipl^^Tipi) / dp2pf'^T{p2) / d6(2 (sin^a)' 



(pip^' +P2Pr' + 2cose2) 't(|pi +P2I) [W{\pi+P2\R)]' B{pi,p2;e2) 



(100) 



where 
B{pi,p2;62 



dpsPi'^ d6l3 sin6l3 (piP3 ^ +P3Pi ^ + 2COS6I3) 
Jo 

dcj>-i [p2P3'^ +P3P2^^ - 2(cose3Cos6l2 + sine3sine2Cos<^3)] ^''^ n(pi,p)n(p2, -p)n(pi,p2) . (101) 
This leaves a 6D numerical integration to be performed. 

• The Chi-squared correction E^2(i?) may be calculated in the same fashion as the contribution Yl^^j^{R); the result is 



15 



(2^) 



^x3/2' 



(4) 



dpipl^^T{pi) / dp2p"'''T(p2) / dSa (sine2)'' 
Jo Jo 

X (pi P2^ +P2Pi' + 2cose2) T(|pi + P2I) [W^dpi + P2|i?)] ' B{pi,p2; 92) n(pi, p)n(p2, -p) , (102) 
where T{p) = r(p)[/3(p) + 21n(l +p/p,„) - 1]"'/^. 

• For the "linearised" LN density field one can perform calculations in a similar way, which lead to equation (jH^), where the 
function (R) = E^"' (R) + J^f (R) is the sum of two contributions originated by the two terms of the linear bispectrum Bg 
in equation (kg): 



.(a),„. 15 A' 



8 (2^)Vo 



dpipi T{piY 



2pIt{p2Y / de2(sin6»2)' 



X (pipa Vp2p/ +2cos6'2) ^ [W"(|pi + P2|i?)]^n(pi,p2 



(103) 



The second contribution is 
^(b)/„N 15 A^ 



8 (2^)*7o 
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dpipl^''T(pif / dp2P^'" / d02(sine2)' 



p. Catelan and T. Theuns 

X {P1P2^ +P2Pl^ +2cOs92y''''' [r(|pi+P2|)]' [Vy(|pi +P2|ii)]'n(pi,p2) . (104) 
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